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Abstract 

Two-dimensional PT-symmetric quantum-mechanical systems with the complex cu- 
bic potential Vu = x 2 + y 2 + igxy 2 and the complex Henon-Heiles potential Vhh = 
x 2 + y 2 + ig (xy 2 — x 3 /3) are investigated. Using numerical and perturbative methods, 
energy spectra are obtained to high levels. Although both potentials respect the VT 
symmetry, the complex energy eigenvalues appear when level crossing happens between 
same parity eigenstates. 
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1 Introduction 



In 1998, Bender and Boettcher introduced the non-Hermitian "PT-symmetric quantum systems 
PQ. This letter ignited a very active field. Hundreds of papers have been published on various 
aspects of this new quantum mechanics. For review, please see Ref. [2] and references therein. 

Most of the studies have been focused on one-dimensional systems, very few have touched 
two or higher dimensions. In 2001, Bender et al studied the complex cubic potential Vu = 
x 2 + y 2 + igxy 2 and Vm = x 2 + y 2 + z 2 + igxyz. Using perturbation theory, they found real 
eigenvalues for the three lowest levels. They used the WKB method to confirm the reality of 
the ground state energy. Their approach also applied to the complex Henon-Heiles potential 
Vhh = x 2 + y 2 + ig [xy 2 — x 3 /3) [3J. In 2002, Nanayakkara and Abayaratne revisited the igxy 2 
interaction with a non-degenerated setup. They computed energy eigenvalues up to n = 4. 
With their choice of parameters, all eigenvalues appear to be real [I]. Later in the same 
year, Nanayakkara obtained the analytic results for igxy interaction and numerical results in 
igxy 3 and igxyz 2 interactions. For all three systems, the author found real spectrum for small 
enough coupling constant and complex one for large coupling. The interaction of igx 2 y 2 is 
also studied [5]. In 2005, Nanayakkara observed avoided level crossings in both the real cubic 
interaction gxy 2 and the complex cubic interaction igxy 2 [6] . The results has been challenged 
in a comment by Bfla et al in Ref. [7]. They pointed out that the quantum system with 
the real cubic potential is ill defined and the complex cubic potential has no avoided level 
crossings. 

In this paper, we compute the eigenvalues of two-dimensional systems to high levels with 
high precision. We study the Hamiltonians associated with the complex cubic potential V12 

H12 = p 2 x + p 2 y + x 2 + y 2 + igxy 2 (1.1) 

and with the complex Henon-Heiles potential Vhh 

#hh = p 2 x +p 2 y + x 2 + y 2 + ig(xy 2 -\x 3 ) . (1.2) 

Both Hamiltonians are not Hermitian. Rather, they respect the VT symmetry and the y 
parity: 

[H,VT]=0, [H,V y ]=0, (1.3) 

where the total parity V, the time reversal T, and the y parity V y are defined as 

V : i —> i; x — > —x; y — > — y; 

T : i -> -i; x->x] y -> y; (1.4) 
V y : x — > x; y -> -y. 

Although the Hamiltonians are complex, the secular equations are real due to the VT sym- 
metry. They can only depend on g 2 . Therefore, both systems are invariant under g — > —g. 

Since the y parity is an unbroken symmetry, we can diagonalize H and V y simultaneously. 
All the eigenfunctions of H are also eigenfunctions of V y with eigenvalues ±1. 

We use different methods to compute the eigenvalues of each system. We found complex 
conjugate pairs when the levels cross between same y-parity states. At higher level, the level 
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crossing happens at smaller coupling constant. This indicates that no matter how small the 
coupling constant is, there are always complex eigenvalues at high enough levels. It shows the 
rich structure in the two-dimensional PT-symmetric quantum systems. One must be very 
careful when generalizing the results in one-dimensional non-Hermitian quantum mechanics 
to higher dimensions. 

The paper is organized as the following: the methods we used are briefly reviewed in Sec. HJ 
the results of different systems are presented in Sec. [3] and Sec. HI we conclude in Sec. El 

2 Methodology 

We use three different methods to study the spectra of Hamiltonians in (11.11) and (jl.2p . 

• The first method is the perturbation theory. Because the energy levels of unperturbed 
Hamiltonian, 

H =pl+p 2 y + x 2 + y\ 

has degeneracy, we must apply the degenerate perturbation theory. After getting the 
eigenvalues as power series of the coupling constant g, we use Pade expansions to extract 
the information from the divergent series. The results for two systems are shown as lines 
in Fig. [1] and Fig. [3j respectively. 

• The second method is numerical computation by the finite-element method (FEM). The 
FEM results for eigenvalues with Vu potential are shown as crosses in Fig. [1] and Fig. [2j 

• The third method is based on the expansion in two-dimensional harmonic oscillator 
(HO) basis [8]. We first analytically compute the non- vanishing matrix elements of the 
full Hamiltonian on the two-dimensional HO eigenfunctions basis. Depending on the 
precision requirement, we truncate the sparse matrix to a finite size. And the eigenvalues 
of the full Hamiltonian can be obtained by numerically diagonalizing this finite matrix. 
This method is a mixture of analytic and numerical techniques. It turns out that this is 
the best way to compute the eigenvalues of the two-dimensional complex Hamiltonians. 
The results are shown as dots in Fig. [1] and Fig. [2] and crosses in Fig. [3J and Fig. HI 

From Fig. [1] we can clearly to see that all three methods provide consistent informations 
about the real part of the eigenvalues. And Fig. [2] shows the two numerical methods are 
consistent on the imaginary part of the eigenvalues. 
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3 Eigenvalues with the complex cubic potential 



Using perturbation theory, we calculate the eigenvalues as power series of the coupling constant 
g. Here are the first few terms for the first 10 eigenvalues: 
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(3.1) 



We actually continue this expansion up to the order of g 40 , then use (20, 20) Pade to extract 
the information from the divergent series. The results are shown as lines in Fig. [TJ From the 
figure, we may see that some Pade expansions have poles and cross each others. 

In the same figure, we show the real parts of eigenvalues from FEM as crosses and from 
HO expansions as dots. The two sets of numerical results clearly agree with each other. And 
surprisingly, up to the crossing points, the Pade expansions fit the numerical results very well. 

Although the Hamiltonian if 12 has VT symmetry, some energy eigenvalues are complex. 
The imaginary parts of the eigenvalues by two numerical methods are shown in Fig. [2j 
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Figure 1: Real parts of the eigenvalues as functions of the coupling constant g with the complex 
cubic potential. Lines are from perturbative expansion using (20, 20) Pade. Crosses are results using 
FEM. Dots are results using the method based on two-dimensional HO basis expansions. Blue disks 
have even y-parity and purple squares have odd y-parity. The system is symmetric under g — > —g, 
only the part with positive g is shown. 
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Figure 2: Imaginary parts of the eigenvalues as functions of the coupling constant g with the 
complex cubic potential. Crosses are results using FEM. Dots are results using the method based on 
two-dimensional HO basis expansions. Blue disks have even y-parity and purple squares have odd 
y-parity. The system is symmetric under g — > —g, only the part with positive g is shown. 



From the Fig. [U we observe two types of level crossings. If the crossing is between different 
y-parity states, the eigenvalues remain real. One example of this type is the crossing between 
the third and fourth lowest states near g — 4. If the crossing is between same ^/-parity states, 
then eigenvalues become complex conjugate pairs. For example, the sixth and seventh lowest 
states cross near g = 1.4. 

We also observed that the level crossing happens in higher levels at smaller value of coupling 
constant. This indicates that no matter how small the coupling constant is, there is always 
complex eigenvalues at high enough levels. There is no critical value, g c > 0, such that for 
\g\ < g c , the entire spectrum is real. 

4 Eigenvalues with the complex Henon-Heiles potential 

In the case of the complex Henon-Heiles potential, a similar patten appears. The energy 
eigenvalues become complex when g ^ 0. There are special features for this interaction. 

• First, there are all rational numbers in every term of eigenvalues in the perturbation 
expansion. 

• Second, there are an accident degeneracy between odd and even ^/-parity. Most of the 
energy levels are degenerated, but interestingly, not all of them. For example, Eio = En, 
E21 = E22, and E30 = E31, but E32 7^ £33 ■ In fact, E32 and £"33 are the same up to the 
order of g 2 , but for the order of g 4 and higher, they separate. This patten remains the 
same for higher exciting states: 

- For even n > 2, E nX = E n2 , E n3 = E n4 , etc. 

- For odd n > 3, E n0 = E nl , E ni = E n5 , etc. 

- For odd n > 3, E n2 7^ E n3 starting from the order of g A . 

• Third, E53 appears to be non-alternating in signs in the expansion. And this is a unique 
case up to n = 9. 
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Here are the first few terms for the first 21 eigenvalues: 
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As in the previous section, we actually continue this expansion up to the order of g 40 , then 
use (20, 20) Pade to extract the information from the divergent series. The results are shown 
as lines in Fig. [3j 

The results from HO expansions are shown as crosses in the same figure. Once again, we 
found the Pade expansions fit the numerical results up to the crossing points. 

The VT symmetry is also broken in this Hamiltonian. Some eigenvalues cross and form 
complex conjugate pairs. The imaginary parts of the eigenvalues by HO expansions are shown 
in Fig. HI 

5 Concluding remarks 

We concluded that both Hamiltonians have complex spectra. The origin of this broken VT 
symmetry is unclear. It could be due to the interaction between Hermitian and non-Hermitian 
potentials as in Ref. [9]. It could also be due to the level contraction instead of level repulsion 
in the non-Hermitian quantum mechanics. 
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Figure 3: Real parts of the eigenvalues as functions of the coupling constant g with the complex 
Henon-Heiles potential. Lines are from perturbative expansion using (20, 20) Pade. Crosses are from 
HO expansions. Even y-parity states are using "+" and odd y-parity states are using "x". 



^*^**ir % ***x, 

++** ** 

++****>< 

: 


*x$** 

*** »*++ »»*»*** 

x*%**+ + -j,!** 4 " 

"4**^ x«»***** ********* - 

: ||i;jjs:"--"' - 


: • 

******* ++**** x 

+ ±**** **** 
- .******* ****** 

; f ****** 


*•*••••» 

*p*j**±+ + , 

'«•=;::*••«;;;;;;;;;••*• ■ 

x****±++ ****»*», 

I* **++ ****»» 

,J S?., ***** 



-4-2 2 



Figure 4: Imaginary parts of the eigenvalues as functions of the coupling constant g with the complex 
Henon-Heiles potential. Only the results from HO expansions are shown. Even y-parity states are 
using "+" and odd y-parity states are using "x". 
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The results indicate that there is no critical value of the coupling constant to have reality 
in this degenerated setup. In the analogue with the igxy interaction studied in Ref. [5], it 
is possible that small enough coupling constant exhibits entire real spectrum in the non- 
degenerated setup. This would be confirmed by a future study. 
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